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O ■ Abstract 

£\j ' We reexamine the production of J ftp and other charmonium states for a variety of 

target-projectile choices at the SPS, in particular for the interesting comparison be- 
tween S+U at 200 GeV/c and Pb+Pb at 158 GeV/c as observed in the experiments 
NA38 and NA50 respectively. The result is a description of the NA38 and NA50 
data in terms of a conventional, quasi-hadronic picture. The apparently anomalous 
suppression found in the most massive Pb+Pb system arises in the present simu- 
lation from three sources: destruction in the initial nucleon-nucleon cascade phase, 
use of coupled channels to exploit the larger breakup in the less bound \ l an d ip' 
states, and comover interaction in the final low energy phase. 
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^ ■ 1 Introduction 



The possible use of J/ip suppression as a signal of unusual behaviour in rel- 
ativists ion collisions, first suggested by Matsui and Satz [1], has attracted 
considerable experimental and theoretical study. Great interest has attached 
to the results obtained by the NA50 collaboration for charmonium produc- 
tion in Pb+Pb collisions at 158 GeV/c: to the early findings presented at the 
Quark Matter 1996 meeting [2] as well as to the startling data later released at 
RHIC'97 [3]. The success of Glauber- like calculations of J ftp production and 
breakup in the p+A and S+U [4-6] systems, coupled with a failure of Glauber 
to provide an equally good description of the apparently accelerated absorp- 
tion in Pb+Pb has been widely interpreted [2,3,5] as a signal of QCD plasma 
creation in these collisions. The very sharp behaviour of the J/ip yield as a 
function of transverse energy E t seen in the later experiment [3] has especially 
attracted attention. 



Preprint submitted to Elsevier Preprint 



9 February 2008 



We attempt to retrace this ground theoretically [7], employing a new, two 
phase cascade approach, described in detail elsewhere [8,9], combined with a 
variation of the Satz-Kharzeev model for production and annihilation of char- 
monium in the initial baryonic collisions. This modeling described below, al- 
lows the coupled-channel aspect of the hidden charm spectroscopy, {ip, x\ V*'} 
to play a more central role. In this first application, we include partons in 
a minimal fashion, to describe for example Drell-Yan production. Hence we 
are testing a 'purely' hadronic description of the anomalous Pb+Pb measure- 
ments. 

It has been pointed out that a hadronic picture might succeed [6] without 
invoking quark-gluon plasma (QGP) creation, if at least part of the seemingly 
anomalous suppression in Pb+Pb could be produced by comover annihila- 
tion, i.e. by interactions of the J /if) with secondary mesons generated in the 
ion-ion collision. The second phase in LUCIFER II, which is a low energy 
cascade, perforce includes the effect of J/ip destruction through such comover 
rescattering. 

In LUCIFER II we have attempted to separate hard and soft processes by 
time scale, see Fig(l), so as to permit partonic and hadronic cascading to 
be joined naturally, in a modular fashion. The separation is effected through 
the use of a short time scale, automatically present at high energies: the time 
Tab taken for the two interacting nuclei A and B to traverse each other in 
the global collision frame. The uncertainty principle allows hard interactions 
involving sufficiently high energy-momentum transfer, i.e. for Q^ 1 < Tab, 
to take place in the first and very rapid cascade. Soft processes involving low 
tranverse momentum are not completed until later. Thus in the initial fast 
cascading the nucleons lose no energy but are still aware of the number and 
nature of the two-particle collisions they have undergone. 

Specifically, the method [9] consists of running the cascade in two stages. The 
first is a high energy fast-time mode in which collision histories are recorded 
and fast processes (Drell-Yan and charmonium production) are allowed to 
occur. Using the entire space-time and energy-momentum history of this stage, 
a reinitialisation of the cascade is performed using elementary hadron-hadron 
data as a strict guide. The final positions and momenta of baryons in the 
first phase, and the number of collisions they suffer are recorded and used to 
generate produced mesons together with their initial momentum and space- 
time coordinates. In the initial ion-ion collision the interacting nucleon paths 
are almost along light-cones. The second cascade begins at Tab-, the time 
of the last nucleon-nucleon collision, with initial conditions specified by the 
reinitialisation, but no secondary interactions are allowed until a formation 
time for produced mesons has passed. The participants in the second phase 
are generic mesons, thought of as of gg-like in character with masses centered 
near M qq = 700 MeV and in the range M q q ~ 0.3 — 1.1 GeV. Generic baryons 
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consisting of qqq are also included and are excited to rather light masses, 
M qqq ~ 0.94 — 2.0 GeV [9]. All the generic hadrons decay via sequential pion 
emission. Normal stable mesons and baryons are also present, and terminate 
the decay chains. See Fig(2) for a pictorial description of the model for basic 
hadron-hadron scattering. 

String-Like Model for Hadron-Hudron Scattering 
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Fig. 1. Time evolution of the ion-ion 
collision. The space-time distribution 
of collisions and decays in hard and soft 
cascades is shown for the minimum bias 
Pb+Pb system. 



Fig. 2. Shown are graphic represen- 
tations of the elements of the model 
for the elementary hadron-hadron colli- 
sion: elastic, single diffractive (SD) and 
non-single diffractive (NSD) scattering. 



Kharzeev and Satz [5] employ a model based on hadronic Glauber theory 
describing production and breakup of the J/ip in ion-ion collisions, to demon- 
strate that such a picture cannot account for the degree of suppression seen in 
Pb+Pb collisions at the SPS. Reasoning similarly, we can make a close com- 
parison of our treatment with their work. The required initial production of a 
cc pair is handled within an effective hadronic formulation both in our work 
and in that of Kharzeev, et.al. . There are, naturally, specific and important 
differences between Glauber theory and a cascade model, and it is partly these 
differences which permit the so-called anomalous suppression in Pb+Pb to be 
explained within a purely hadronic framework. 

The overall degree of suppression in Pb+Pb, insofar as it differs from earlier 
work [5,6], results from a combination of effects; these are baryonic, coupled 
channel and comover in kind, with substantial contributions arising from both 
phases of the cascade. There are potential unknown variables: the production 
and dynamic time evolution of each charmonium state, the breakup proba- 
bilities against both baryons and mesons, the density of secondary mesons. 
This last is to a large extent predicted by the cascade, which must agree with 
actual inclusive final state meson and baryon distributions. There also exist 
constraints on the basic charmonium variables. The production is in principle 
determined in elementary nucleon-nucleon collisions, the baryonic breakup in 
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p+A collisions. The ipn breakup cross- sect ions are not known directly from 
any measurements. 

We leave in this conference proceeding report the details of the cascade archi- 
tecture to the perusal of a previous publication [9] . 



2 Coupled Channel Model for Charmonium 

The treatment of the hidden charm cc mesons within a 'purely' hadronic 
code presents some problems, perhaps not fully solvable within the effective 
hadronic treatment of such states. We do not deviate much in spirit from the 
work of previous researchers [4-6], but the devil lies sufficiently in the details 
to produce some quantitative effects. The production of charmonium mesons 
is almost completely limited to that coming from nucleon-nucleon collisions 
at the highest energies, i.e. in the initial high energy cascade, not by fiat but 
by the greatly reduced collision energies in the second phase. Destruction of 
the charm meson precursors, in contrast, can occur in the first baryonic phase 
and also later in collisions with generic mesons and baryons in the second, 
low energy phase, i. e. on comovers. It is in the destruction of the charmonium 
states we differ most, ascribing a more direct role to the presence of the higher 
mass x an d ip' mesons, for which in fact breakup is far easier. We include in 
Fig(3) a level diagram showing the relevant charmonium states to make the 
picture as clear as possible. 

In the actual calculations of the above cited references both production and 
breakup are treated as instantaneous. There is no J/ip formation time in 
the high energy phase or its near cousin, the Glauber or eikonal modeling. 
Kharzeev et al. [5] in fact justify such a choice by referring to microscopic 
production of charmed quark pairs, the subsequent formation of a preresonant- 
resonant state from which all charmonium mesons emanate, and the breakup 
in relatively hard scattering by gluons radiating from nearby nucleons. 

We imagine that the primordial cc pairs are originally produced essentially in 
plane wave states. Clearly, both singlet and octet color states are involved. 
This view is reasonable given the predominance of open charm production 
over hidden charm production in free space NN collisions. We further suppose 
that in elementary NN collisions the cc pair eventually coalesces, with a state 
dependent probability, into a J/ip, ip' or x- The time which elapses will be 
determined by the size of the bound state and the probability that a transition 
occurs. The probability of formation will depend critically on the relative 
momentum of the coalescing pair as well as on their spatial separation. 

Therefore, whether one sees the early evolution of the eventual charmonium 
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Fig. 3. Charmonium spectroscopy in- 
cluding higher mass states which are 
significantly produced in pp and which 
feed strongly to the Jip. Electromag- 
netic and hadronic decays of \ l ( a 
weighted average) and ip' are both in- 
cluded in the indicated branching ra- 
tios. The production ratios are sug- 
gested by direct measurement. 



Fig. 4. Production of J ftp from pp as 
a function of energy. The data appear 
as points and the fit as a line. The irp 
cross-section is also known, and in fact 
is very similar to that for pp, but rarely 
plays a role with production generally 
significant only at the highest energies. 



as a preresonant state or as plane waves may be immaterial. Given the small 
size of the J/0 as opposed to the much larger ip' and Xi the separation of the 
c and c in the plane wave picture could equally well serve as a distinguishing 
feature. What can differentiate our calculation from the earlier models is the 
possibility of transitions between charmonium states dynamically in a nuclear 
environment. Certainly, the \ l ) X 2 an d X 3 states are produced considerably 
more copiously in basic pp collisions, with perhaps as high a ratio as xjip = 
4 — 5 [17], and they decay appreciably into J/ip, with branching ratios in the 
range T^/V ~ 12 — 25% [18]. The ip' also feed some 57% into J/ip. It follows 
that one cannot ignore their presence. 

This point becomes even more significant when one considers what the breakup 
probabilities for the higher mass charmonium states are likely to be, either 
in the fast or slow cascades. These heavier objects are considerably larger 
spatially and might well have total cross-sections on baryons or mesons pro- 
portional to the square of their colour dipole radius [19]. One of our conclu- 
sions will be that a considerable portion of the anomalous suppression seen in 
Pb+Pb, even for quite large impact parameters, is a result of breakup in the 
higher charmonium states and an extinction of their free space feeding down. 

What extra parameters has our model introduced relative to other treatments? 
We introduce a reaction matrix for charmonium states Rij which permits 
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transitions between the states as well as diagonal, breakup, elements. Unitar- 
ity constrains this matrix. The diagonal elements should be present for other 
practitioners also but are not, in general, since only one preresonant char- 
monium state is usually considered. The quantitative reaction matrix is 
specified in ensuing sections, but is in this work limited to diagonal. 

The formation time for secondary mesons, t/, controls the density of comovers 
at the onset of charmonium destruction. A reasonable choice for this parameter 
is tj ~ 0.5 — 1 fm/c, in fact consistent with the production of n mesons at the 
SPS [9,20,21]. The effective formation time is actually somewhat longer, since 
it is increased by the duration of the fast cascade, i.e. r e ff = ry + Tab/2- 

The energy dependence of the elementary J ftp production cross-sections is 
shown in Fig(4). The sharp dependence of aj/^ on energy near the SPS values 
= 17 — 20 GeV implies that virtually all production occurs in the high 
energy phase. 



3 Drell-Yan 

The high energy phase, designed to record the initial interactions of the nucle- 
ons in the two colliding nuclei also provides the basis for our estimate of mas- 
sive dilepton production, i.e. Drell-Yan, an important side of the quandary 
we faced at the start. We limit ourselves to the canonical FNAL [11] p+A 
measurement at 800 GeV/c, but in fact the method of calculation guaran- 
tees agreement with the lower energy p+A and A+B collected by NA50[2]. 
Drell-Yan is generally considered to be calculable perturbatively for dilepton 
pairs with masses in excess of M MM = 4 GeV. Production in the short time 
defined by such masses proceeds without energy loss and leads to the linear 
A-dependence shown in Reference [9]. 

Any cascade which does not correctly describe the A dependence of Drell- 
Yan is in danger of producing spurious charmonium suppression by means of 
premature energy loss. This point cannot be emphasized enough. 



4 Charmonium Suppression in Nuclear Collisions 

4-1 Minimum Bias: p+A and comparison to Glauber 

We begin with the suppression in p+A for which comovers play little role. 
Even here, the first stage high energy cascade does not suffice for an accurate 
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description, some of the suppression on baryons occurs only in the second 
stage, as slow J/^'s emerging from the interaction region are caught by nu- 
cleons, or interact at low energy in the target. To facilitate a comparison with 
the cascade we have made our own calculations with the Glauber formalism 
[7,5]. 

It is instructive to extend this comparison to A+A collisions to demonstrate 
that even Glauber does not reproduce the canonical power law, implied in the 
experimental descriptions [2,13,20], supposedly arising from purely baryonic 
breakup. These results are displayed in Fig(6) for the J/ip. The absorption 
cross-section taken to reproduce the p+A observations at 800 GeV/c [11], 
Cabs ~ 7.0 mb , is equally successful for the lower SPS energies. 

A second comparison can be made using the coupled channel modeling, which 
we employ in the full LUCIFER calculations below. The relative production 
of the different charmonium states is taken so as to reproduce the pp data 
from the ISR [17] for the x to J /if) ratio, i.e. xl^ ~ 4.5, and for appropriate 
ip' production [22]. Our final results are rather insensitive to a choice in this 
range, since a decrease could easily be compensated by a small transition 
matrix element between J/ip and x- 

A lesson, key to our development, is that the coupled channel model repro- 
duces the Glauber result for J/ip, using a smaller direct breakup cross-section, 
Cabs(J/ip) = 5 — 6.0 mb, but including indirect destruction via the considerably 
larger <J a bs{x) = 3a a f, s (J/^) for x an d perhaps higher for if)'. The increased 
spatial sizes of the higher states strongly support the use of larger absorption 
cross-sections. One observes that pure Glauber theory and the first stage nu- 
cleon cascade, both produce deviations downward from any power law, and 
thus gathers there is a little bit of 'anomalous' suppression even in a bare 
bones, no comover, theory. 

4-2 Suppression in A+B Collisions 

To complete the picture one must allow the soft cascade to go forward for 
ion collisions where the production of mesons becomes significant. There are 
two sets of data to be considered: first, minimum bias J/ip cross-sections as a 
function of the product Ax B of nuclear atomic numbers, and second the ratio 
of J/ip yield to Drell-Yan yield as a function of centrality, or more specifically 
transverse energy E t . 

Our results for minimum bias are displayed for the combined effect of both 
cascade phases in Fig(5). The anomalous suppression in Pb+Pb is well repro- 
duced by the totality of our two step, but otherwise conventional, hadronic 
dynamics. Part of the additional suppression in Pb+Pb relative to S+U al- 
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ready arises from the high energy cascade, coming from the increased x an d 
ip' breakup in the more massive nuclear collision. But a considerable differen- 
tial suppression arises from comovers, both mesonic and hadronic, some 40% 
of the difference between S+U and Pb+Pb. Part of the anomaly however, is 
perhaps illusory in view of the 'curving down' seen for large in Fig(6). 

JA\l Suppression in pA and AB 
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Fig. 5. The whole range of yields for 
J/ip from pp and p+D to Pb+Pb cal- 
culated in the cascade and compared 
to SPS measurements at various en- 
ergies. The absolute theoretical values 
are obtained by normalisation to nu- 
cleon-nucleon. 



Fig. 6. Comparison for A+A between 
Glauber and cascade, the latter in a 
purely J/ip mode and both calculations 
employ a^r = 7.3 mb. The deviation 
from a power law is apparent for large 
Ax A. A hard sphere form is used for 
the nuclear density. 



The calculated minimum bias ip' suppression is compared to data in Fig (7). 
The strong drop occasioned by the large increase from p+W to S+U or Pb+Pb 
is clearly present in the theory. As is evident in this figure the ip' breakup 
strength inferred from p+A proves sufficient for both S+U and Pb+Pb. 



The breakup cross-sections in these simulations are 6.6, 20.0 and 25.2 mb 
for the ip, \ % and ip' respectively. These represent absorption in charmonium- 
baryon collisions, and are reduced by the constituent quark factor 2/3 in ip- 
meson. Variation of these meson-meson cross-sections upwards to full equality 
with charmonium-baryon leads to ~ 0.5% change in the overall J/ip suppres- 
sion for Pb+Pb (see Fig(8)). 

This perhaps surprising non-linearity for the J/ip interaction arises because in 
our model this state is intrinsically tied up with the higher states. Introducing 
off diagonal elements Rij would produce a family of solutions. We have left well 
enough alone; the present few modeling parameters, now mostly determined 
independently of S+U or Pb+Pb data, surely having produced already an 
adequate description of observation, and any attempt to force a statistically 
accurate theory of the NA50 data, though possible, is not justified in such 
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Fig. 7. Comparison of experiment vs 
simulation for ip'. The Pb+Pb data 
from NA50 was rescaled to 200 GeV/c 
by the collaboration. The S+U data 
is taken from NA38. The cascade 
calculations, again normalised to nu- 
cleon-nucleon, reproduce the observed 
behaviour for p+A and the sharp drop 
in the ip' to J/ifi branching ratios for 
the massive nuclear collisions. 

schematic models. 



Fig. 8. Varia- 

tion of J/ip suppression with the char- 
monium-meson cross-sections. We use 
2/3 as the ratio to charmonium-baryon 
for the calculations in the paper. The 
small variation with cross-section size 
is surprising, and indicates a saturation 
is taking place in the charmonium sup- 
pression. 



The survival probabilities for the J/ip in S+U are 0.50 and 0.67 in the hard 
and soft cascades respectively. The same figures for Pb+Pb are 0.42 and 0.775. 



4-3 Centrality: Dependence on Transverse Energy 



Perhaps the most striking features of the NA50 [3] measurements are contained 
in their plot of J/ip suppression vs E t . Unlike the existing Glauber calculations 
of transverse energy the cascade provides a built in E t scale, which does not 
necessarily agree exactly with the experimental determination. NA50 plots for 
J /if) and Drell-Yan show E t from neutral energy within the pseudo-rapidity 
range rj = 1.1-2.3. To establish a calibration from LUCIFER II we first refer 
to their earlier Pb+Pb results [20] using a more central rapidity range r) — 2.1- 
3.4, and including both electromagnetic (neutral) and hadronic calorimeters to 
estimate E t . This comparison is shown in Fig(9) [20], and indicates that LU- 
CIFER II, with standard parameters [9], provides a reasonable representation 
of the measurements. The small discrepancy between cascade and experimen- 
tal endpoints, some 10 — 15%, should be kept in mind when examining the 
NA50 charmonium data. 
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Figures (10) and (11) display the results of simulations for the two massive 
ion-ion collisions. The magnitudes use the calculated survival probabilities, 
normalised by the pp or p + D experiments. The rather low E t value at which 
the measured J /if) suppression becomes pronounced more or less obtains in 
the simulation, and the low level of J/0's for higher E t is reproduced. The 
theoretical errors shown are only indicative of the monte-carlo simulation and 
given the normalisation to Drell-Yan ratios for nucleon-nucleon should include 
some further systematic normalisation error. Appropriate integrals of the E t 
plots are, however, consistent with the minimum bias results in Figure (5) for 
both S + U and Pb + Pb. 

The results reinforce the perception already created by the comparison with 
the minimum bias data. The hadronic two-step cascade is capable of describ- 
ing the charmonium yields for J /if; and if;' as well. The beginning of strong 
suppression in J /if; at close to peripheral collisions, i. e. at low E t , is a reflec- 
tion of the role the heavier charmonium states play. The scale used for central 
S + U is just that obtained from the cascade. For Pb + Pb, where NA50 used a 
more peripheral cut to obtain the neutral energy the theoretical cutoff is close 
to 125 GeV, somewhat lower than say that quoted by the experiment [2,3], 
and a scaling factor of ~ 1.25 has been applied to the theoretical transverse 
energy, justified by subsquent discussions with the experimentalists [23]. The 
original theoretical scale is, however, closer to that inferred from the NA49 
measurements Reference [20]. Of this factor of 1.25, some 10 — 15% is, as we 
indicated, attributable to differences in the NA49 E t calculated and measured 
scales. 



5 Conclusions 

It appears that a conventional hadronic explanation of the minimum bias and 
central J /if; and if>' suppression in A+B collisions is possible. This is accom- 
plished here with a cascade, not specially tuned for the charmonium sector 
alone, but consistent with soft energy loss processes, both meson and proton 
spectra, and Drell-Yan data. There are parameters in the model, notably the 
meson formation time, the generic meson decay constant, and certainly the 
elements of the charmonium reaction matrix. However, we have not made use 
of all of this freedom in obtaining the main results. Indeed, the only parameter 
not constrained by independent information is the charmonium breakup cross- 
section on mesons. The variation of suppression probability for a reasonable 
range of this variable, as seen in Figure (8), is very small. 

Some authors [5] suggest theoretically that the J /if} total cross-section on 
mesons must be drastically smaller than the ~ 4 mb we use for breakup, but 
this argument is disputed by others [19], and it seems unlikely. In any case 
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Fig. 9. Transverse energy distributions 
from LUCIFER compared to exper- 
iment, for NA49. A broader range 
of pseudo-rapidity has been used and 
both neutral and charged energies are 
present. The calulation has underesti- 
mated the measured endpoint by some 
10-15%, which must be taken into ac- 
count in considering charm suppres- 
sion. 



Fig. 10. Comparison between the cas- 
cade and NA38 transverse energy de- 
pendence for J ftp, for S+U. One notes 
that the pseudo rapidity range used for 
the neutral Et in the case of S + U is 
broader, i] = 1.7-4.1, than in NA50. 
The calculated E t spectrum shown was 
obtained with no scale factor. 



breakup on mesons occurs at quite high relative energies where the very small 
cross-sections suggested by these authors, even if valid, would surely not ob- 
tain. Our direct experimental knowledge of the total and partial cross-sections, 
including any energy dependence, for J/ip or other charmonium mesons on the 
lower mass mesons is of course very limited. 

What then has been learned about excited, dense, nuclear matter from the 
reduction in J/^'s? Our earlier calculations [9] for a broader range of processes, 
suggested that very high baryonic and mesonic energy densities were achieved 
in central Pb+Pb interactions, ps ~ 4-5 GeV/(fm) 3 and p# ~ 3 GeV/(fm) 3 
respectively and that these densities persist for quite long times r ~ 3-5 fm/c, 
in the cm. frame. Thus appropriate conditions for possible 'plasma' creation 
exist in the most massive collision. This matter density has been sensed in 
the theoretical comover breakup of charmonium, both J/ip and ■?/>', more so 
for Pb+Pb than for S+U. But should our model stand the test of time, and it 
has explained a good portion of existing data at the SPS, then the case for a 
non-conventional explanation is hard to establish as yet. It is always of course 
possible that partons actually are playing a less passive role than portrayed 
by hadronic modeling and that especially high gluon densities are achieved in 
the initial phase. Our simulations do not rule out the creation of some form of 
partonic matter in ion collisions at SPS energies. They only make the necessity 
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Fig. 11. Comparison between the cas- 
cade and NA50 transverse energy de- 
pendence for J/V>. There are no discon- 
tinuities, of course, in the LUCIFER 
yields, but the general shape is repro- 
duced. See the text for a discussion of 
the E t scale. 

thereof less compelling. A microscopic, non-hadronic treatment of the internal 
strusture of the charmonium states might alter the entire picture, and may be 
necessary at higher energies. 
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